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SUMMARY 


_The aim of this investigation is to estimate the distribution function of lattice vibrations at 


0°K in KCl. Earlier Iona (1941) determined it from the elastic constants and the residual ray 
frequency measured at room-temperature using the classical Born model. This computation 
is recalculated with new empirical constants measured at 293°K. Since the residual ray frequency 
is not known at 0°K, the Born model according to Iona cannot be used in the calculation of the 
distribution function at 0°K. Therefore we use a modified model where the classical ionic charge 
‘is replaced in the Coulomb interaction by an effective charge, which Lundqvist (1955) studied 
in a quantum mechanical investigation. The specific heat c, and the Debye temperature 6, are 
computed from this distribution function. The agreement between the calculated values of 0p 
and the experimental values at low temperatures is not so good as that obtained by Iona. In a 
preliminary calculation the agreement seems to be better if the polarization, here excluded, is also 


taken into account. 


Introduction 


The distribution function of thermal vibrations in ionic crystals has earlier been 
calculated in connexion with investigations of the specific heat at low temperatures. 
Kellermann [10] estimated it for NaCl and Iona [8] for KCl, both using the classical 
Born model. More recently Lundqvist [12] investigated the possibilities of calculating 
the vibrations from the quantum mechanical expression for the binding energy in 
the Heitler-London approximation. 

Some inaccuracy seems, however, to exist in Iona’s and Kellermann’s calculations, 
because they used the elastic constants measured at room-temperature instead of 
those measured at 0°K, which have to be used in order to compare the theoretical 
results of the specific heat c, with the experimental ones. The empirical constants 
have later been more carefully measured. On account of this a closer investigation 
of their results seems to be necessary, because their calculations have often been 
referred to in experimental measurements of c, at low temperatures. 

Our investigation follows in the main the preceding ones. The forces acting upon 
the ions are regarded as central forces and the KCl lattice as a simple cubic lattice. 
The purpose is to recalculate Iona’s results with better empirical constants. The 
coefficients in the equation of motion are determined by the elastic constants, the 
equilibrium condition and the residual ray frequency, by taking into account the 
repulsive forces from the nearest and the next nearest neighbours. 
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0 | can be expected in the region where the temperature is below 
it is, however, impossible to use the Born model according to Iona, 


the residual ray frequency has only been measured at room-temperature. In o1 

calculate a distribution function at 0°K we use a simplified approach by introdu 

an effective charge e* and consider the repulsive forces only between the neares 

4 neighbours. — is } nee Asie o> 

; The effective charge has been investigated by Szigeti [18] and more recently by 
Lundqvist [12] and [13]. Although the concept of an effective charge was introduced 

only in connexion with the long polar oscillations, it follows from Lundqvist’s 
investigation that to a certain degree of approximation the same effective charge 
can be used for all modes of vibrations in the special case of KCl. 

Computations by this method are first done-at room-temperature in order to 
make a comparison with the calculation according to Iona’s scheme. The comparison 
shows that the transverse branches change very little and the longitudinal branch is 

S cut down to a better agreement with the value calculated from electromagnetic 
theory. 

Lastly we estimate the distribution function at 0°K with the method sketched. 

From the distribution function c, the Debye temperature 0, is determined. 6, is not 

in good agreement with the experimental result when the temperature is above © 

10°K, which may be due to the polarization which is not taken into account here. : 

A preliminary calculation shows that the frequencies seem to decrease generally | 

if it is included. Thus the distribution function is displaced to smaller frequencies, — 

the specific heat increases, and the Debye temperature decreases, in better agreement — 

with the experiment. 


The frequency spectrum 


The KCl lattice is assumed to be a simple cubic lattice with equal masses for K* 
and Cl-. The forces acting upon the ions are regarded as central forces. In this case 
the Cauchy relation c,, =c,, is fulfilled. For KCl and NaCl the deviation is not so 
large. Lowdin [15] first showed that this deviation may depend on many-body 
forces which he separated and found to be of importance in explaining the difference 
between cy. and c,4. In our calculation the mean value of c,. and c,, is used. 

The potential energy is expanded in powers of the displacerhents from the equilib- 
rium position and the second order term is put in the equation of motion. Looking 
for plane wave solutions with the frequency w and the amplitude U we obtain 


mo? U,=—> 4.3 U, (three equations), (1) 
¥ 


where m is the mass and the A,, are given from the derivatives of the potential 
energy of the ion in the direction of # and y. 

We divide the potential energy into two parts, one electrostatic and one due to the 
repulsive forces. Accordingly the coefficients A,, are divided into two parts A,, = 
Aty + Ajy. Iona has calculated A’, from a formula given by Waller for a simple 
cubic crystal. A}, can be estimated from the force constants P and Q according to 
Born [1] and Born and Géppert-Mayer [2], when the elastic constants are known. 
One relation is obtained from the equilibrium condition. When the repulsive forces 
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ations an ive ck * is introduced. It has been in ted 
18] and more recently by Lundqvist [12, 13]. Szigeti deduced-an expres- 
‘it using macroscopic theory and obtained a formula containing empirical 
nts. For the alkali halides he found that e*/e is smaller than unity and for 
KCl equal to 0.80 when the empirical constants at room-temperature are used. _ 
~ Lundqvist [12] considered the quantum mechanical expression for the binding 
energy in the Heitler-London approximation. He used the formulas for the binding 
energy of an ionic crystal, which are based on a work by Hylleraas [7] and developed 
Landshoff [11] and later discussed in great detail by Lowdin [15]. In the approxi- ~ 
mation considered by Lundqvist, the potential energy contains both two-body 
‘interaction and also many-body interactions. He pointed out that these many-body | 
terms lead to the introduction of an effective charge. He derived a formula, which is 
related to the elastic constants c,, and ¢,,, from an approximate expression for the 
three-body potential and obtained the value 0.86 for KCl when the elastic constants 
at 0°K are used. Later [13] he improved the approximation for the three-body 
_ potential, giving rise to certain correction terms in the effective charge. The expres- 
sions for the coefficients were derived only for the long polar vibrations but from 
_ these expressions the behaviour at other frequencies is clear. If the residual ray 
_ frequency were known at 0°K it would be possible to deduce a good value for e* 
from it, but so far it has only been measured at room-temperature. When the Cauchy 
relation is valid one gets a good result for e* from c,5, but the larger the deviation, 
_ the worse is the resuit. Although the effective charge generally is a rather complicated 
effect, one can take it into account for the special case of KCl in an approximate 
" way simply by replacing the classical ionic charge by an effective charge in the Cou- 
_ lomb interaction, because the deviation from the Cauchy relation is small for KCl. 
eFor KCl at 0°K ¢c,,=0.6 x 10% dyn/em? and c,, = 0.669 x 10" dyn/cm?. This . 
deviation may depend on many-body forces, which are neglected in our calcula- 
tions, and the mean value of c,, and c,, is used. The effective charge is estimated 
from the expression for c,. when the repulsive forces from only the nearest neigh- 
bours are taken into account. In that case only the electrostatic part of the coefficients 
contributes to c,, and the expression for it is according to Léwdin [14] 


*2, 
3¢ 0.2318. (2) 
a 


C12 = 


The value obtained for e*/e is 0.86 at 0°K and the same at 293°K." 
In one of our investigations we use Szigeti’s value and consider the interaction 
between the nearest and the next nearest neighbours. In this case we have no theo- 
retical justification for the calculation and we wish only to see how the frequencies 
are affected. 
The calculations are done for five different cases: 


1 This formula is approximately correct only for a static lattice. For the present work, however, 

it will be assumed to be valid for arbitrary temperatures, in order to give the correct elastic 

_ waves. Corresponding assumptions seem to have been made in all calculations of lattice vibrations. 
As far as we know, no theoretical justification for such assumptions has been given. 
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-, TENERZ, Calculation of lattice vibrations in KCl 


A. At 293°K with the repulsive forces from the nearest and the next nearest 

neighbours and e*/e = 1. ; 
B. At 293°K with the repulsive forces from the nearest neighbours and e* /e = 0.86 
C. At 0°K with the repulsive forces from the nearest neighbours and e*/e = 0.86. 


D. At 0°K with the repulsive forces from the nearest neighbours and e*/e=ia 


E. At 293°K with the repulsive forces from the nearest and the next nearest 


neighbours and e*/e = 0.80. 
See Table 1. 
Table 1. Values of the frequencies. 
) PRES SR an ee SPRUE Be 
A B Cin, D 
oO oO oO oO oO oo oO od o o oO o 
pect apts: dhe th eS le Behe eae eae eee 
e re re = || Comal = Lom) ee = _ Loom! 
x x x x x x x x x x x x 
Sega’ Sees aE eet rg rg rg Heme ty aba years. 
| 
0] 0| 4 || 1.09 | 0.44] 0.44 | 1.09 | 0.48 | 0.43 | 1.21] 0.43 | 0.43 | 1.21} 0.51 | 0.51 
0| 0 | || 1.99 | 0.83 | 0.83 | 2.00 | 0.80] 0.80 | 2.22 | 0.81 | 0.81 | 2.22 | 0.95 | 0.95 
0| 0 | 2 /\2.59] 1.10 | 1.10]| 2.60} 1.06 | 1.06 | 2.89 | 1.08 | 1.08 || 2.88 | 1.25 | 1.25 
0| 0] 1|/2.79| 1.20] 1.20) 2.81 | 1.16] 1.16 | 3.12] 1.17] 1.17 | 3.11] 1.36] 1.36 
0] 4] 4 1.30] 1.01| 0.68, 1.29} 1.01 | 0.65 | 1.39] 1.14] 0.66 | 1.45] 1.10] 0.77 
0| +| 4 2.04] 1.30] 1.07 | 2.04 | 1.28 | 1.02 | 2.25 | 1.39 | 1.03 || 2.27 | 1.44] 1.20 
0| | 3 | 2.59] 1.49 | 1.38 | 2.60] 1.46 | 1.30) 2.89 | 1.56 | 1.32 |) 2.88 | 1.66 | 1.53 
0| 4| 1 | 2.78} 1.57] 1.50/) 2.79 | 1.53 | 1.41 | 3.11 | 1.62 | 1.42 | 3.09 | 1.75 | 1.66 
0| 4] 4 | 2.26] 1.89 | 1.55 | 2.23 | 1.88] 1.43 | 2.43] 2.11 | 1.45 || 2.50 | 2.07 | 1.68 
0] 4] | 2.62 | 2.11] 1.94 2.61 | 2.09 | 1.79 | 2.89 | 2.31 | 1.80// 2.89 | 2.32 | 2.10 
0| 4| 1 | 2.76] 2.18} 2.10/| 2.76 | 2.14] 1.92 3.07 | 2.35 | 1.94 || 3.04 | 2.39 | 2.26 
0| 2] & | 2.67] 2.50 | 2.44 | 2.62 | 2.51} 2.21 | 2.93 | 2.78 | 2.24 || 2.93 | 2.74 | 2.60 
0| #| 1 | 2.73 | 2.63 | 2.59 || 2.72 | 2.57 | 2.38 | 3.04 | 2.85 | 2.41 || 2.99 | 2.83 | 2.80 
0} 1| 1) 2.84 | 2.72 | 2.72 || 2.56 | 2.71 | 2.71 | 2.59 | 3.02 | 3.02 || 3.01 | 2.97 | 2.97 
3] 4| 41.57] 1.07] 1.07// 1.53 | 1.07] 1.07] 1.63 | 1.19 | 1.19 || 1.75 | 1.19 | 1.19 
4/4] 4] 2.18] 1.57 | 1.23] 2.14] 1.52 | 1.22 || 2.83 | 1.64] 1.33 || 2.41 | 1.71 | 1.37 
+| 4| $2.63 | 1.96 | 1.38 || 2.62] 1.85 | 1.36 | 2.90] 1.95 | 1.46 || 2.91 | 2.13 | 1.54 
4/4] 1] 2.77] 2.12] 1.44]! 2.78 | 1.99 | 1.42 | 3.09] 2.07 | 1.52 | 3.07 | 2.30 | 1.61 
1) 4 4 | 2.57] 1.90| 1.62 || 2.47 | 1.89 | 1.57 || 2.64 | 2.12 | 1.69 || 2.81 | 2.08 | 1.77 
4/4 | 2 2.85 | 2.36] 1.81) 2.71] 2.29] 1.76 2.99 | 2.50| 1.89 | 3.09 | 2.57] 1.99 
+| 4| 1 || 2.76 | 2.75 | 1.86 || 2.74 | 2.57 | 1.82 || 3.06 | 2.71 | 1.97|| 3.02 | 2.95 | 2.04 
| | 2 |/3.10| 2.49 | 2.33 | 2.91 | 2.48 | 2.24//3.13 | 2.78 | 2.38 | 3.32 | 2.73 | 2.53) 
+| 2] 1// 3.13] 2.72 | 2.41 || 2.89 | 2.71 | 2.35 || 3.07 | 3.03 | 2.53 || 3.32 | 2.98 | 2.62 || 
bf} 1) 1 a 2.71 | 2.71 | 2.87 | 2.69 | 2.69 | 2.95 | 3.01 | 3.01 | 3.35 | 2.96 | 2.96 
4] 4] [3.11] 1.93 | 1.93 | 2.92 | 1.92 | 1.92// 3.09 | 2.15 | 2.15 || 3.34] 2.11 | 2.11 | 
k| 4] 2|/3.43 | 2.42] 1.96) 3.18 | 2.40] 1.95 | 3.37 | 2.67 | 2.17]| 3.64 | 2.64 | 2.15 || 
! 3 ee 2.72 | 1.98 | 3.16 | 2.71} 1.96 | 3.33 | 3.03 2.18 | 3.63 | 2.98] 2.16. 
$ | 2|/3.80| 2.48 | 2.31] 3.48 | 2.46 | 2.28 3.68 | 2.76 | 2.53 | 3.99 | 2.70 | 2.51 
4/2] 1) 3.86] 2.70] 2.37] 3.50] 2.68 | 2.35 | 3.68 | 3.00 | 2.62 | 4.02 | 2.94 | 2.58) 
b] 1) 1 3.90] 2.67 2.67 | 3.50 | 2.66 | 2.66 3.67 | 2.98 | 2.98 | 4.05 | 2.92 | 2.99 
$/2| 2 | 4.28] 2.48 | 2.48 | 3.90] 2.44] 2.44 | 4.12 | 2.74| 2.74 || 4.47 | 2.68 | 2.68 | 
£) 2] 1) 4.39} 2.68 | 2.47|/ 3.97 | 2.64 | 2.44) 4.19| 2.97| 2.73 | 4.55 | 2.90 | 2.67 | 
$| 1] 1 | 4.50] 2.67 | 2.67/| 4.04 | 2.63 | 2.63 4.26 | 2.96 | 2.96 | 4.65 | 2.88 | 2.88 
1] 1] 1|/4.72 | 2.67 | 2.67 || 4.25 | 2.62 | 2.62|/ 4.49 | 2.95 2.95 | 4.87 | 2.87 | 2.87 | 


Ha= Pal. (x=, y, 2,) 


280 


at 293°K. 


a ic heat the foll wing relation is valid: 
ae : — 

oe Ae a ty NB 
q rons (7o) eur Ln 
Ee” ,  O=kNf we — to) dy, (3) 
(e*T — 1) 


where N is the number of the oscillators, & Boltzmann’s constant, 4 Planck’s con- 
stant, v are the frequencies of the oscillators, and f(y) is the distribution function of 


_ lattice vibrations. 


order to use this formula it is necessary to know the distribution function, 


_ which is computed in the following way. The frequencies are determined from the 


equations (1). Because the KCl lattice is regarded as a simple cubic lattice, the 


_ determinant gives a third degree equation in w? from which w can easily be found. 
: From the determinant three positive, in general different, values of w are found. In 
_ Table 1 these are arranged in three branches of vibrations, one containing the highest, 
_ one the intermediate, and one the lowest frequencies. In order to calculate the 
_ distribution function we interpolate the frequency midway between two calculated 
- points. This is possible by constructing three smooth surfaces, one for each branch 


+ 


4 


in the ~, my, pz w space, where 7, P, y, are the components of the wave vector. Thus 
we obtain the frequencies at the points where y,=nz/8, p, =m/8, p, =12/8 


(n,m, l integers). 


F 


The whole frequency range is divided into 24 intervals, and the number of cal- 
culated frequencies in each interval is counted. Finally we normalize the distribution 


~ function so that 


[f0)dv=3. (4) 


- Table 1 gives the frequencies in the different cases. 


The frequencies are calculated at 293°K with the Born model according to Iona. 
A comparison between Iona’s results and ours shows how a change of the elastic 


~ constants alters the result. Iona used the elastic constants measured by Férsterling 


[6] which are c,, =3.767 x 101! dyn/cm? and c,, = 0.633 x 104% dyn/cm?, while we 


use ¢, = 4.00 x 10 dyn/cm? and c,, = 0.616 x 104 dyn/cm?. Our frequencies are 


in general somewhat higher than Iona’s. 
The aim of the investigation is to determine the distribution function at 0°K, 
but as long as the residual ray frequency is not known at 0°K, we cannot use the 


Born model according to Iona. For this reason an effective charge is introduced 


which is estimated from c,, and the repulsive forces between nearest neighbours only 
are taken into account. A comparison between this calculation at room-temperature 
and the preceding one shows that the agreement for the transverse branches is good 
and that the limiting frequency for the longitudinal branch is in better agreement 
with the value 3.96 x 1013, which is obtained from the dielectric constants. It con- 
firms the assertion that the frequencies will be nearly the same when only the nearest 
neighbours are considered for the repulsive parts of the coefficients and when the 
effective charge is computed from the experimental values of ¢19. For this reason a 
calculation at 0°K is done with an effective charge and the repulsive forces from the 
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Fig. 2. The distribution function 
calculated from the frequencies of 
table 1 C. 


0 1 2 3 4 5 6 7x1012 sec-3 


nearest neighbours only. We see that the limiting frequencies at 0°K are in good 
agreement with the values obtained by Lundqvist. 

Figs. 1, 2 and 3 show the distribution function in three different cases. From them 
one can see that the highest maximum of the functions moves to higher frequencies 
when the temperature decreases. The maximum lies at the residual ray frequency, 
and when the temperature decreases, the residual ray frequency moves to higher ' 


282 


of ote 


— Tie th) 


. ta %g 
» . ‘ 
’ : 4 
a 
\ : 
ei oon 
3 


5 6 71012 sec-1 


Fig. 4. Debye’s characteristic temperature. 


so - + —Clusius, Goldmann and. Perlick [4]. 
: - — Keesom and Clark [9]. 
I — calculated from fig. 2. 
II — calculated by Iona. 


_ values, and so the peak too is moved.! We can also see that the first peak is not as 
_ high at 0°K as at 293°K. This peak has a great importance in the calculation of c, 
at low temperatures. 

A comparison between the zero-point energy calculated with Debye’s distribution 
function and that calculated with ours at 0°K shows that the Debye model gives 
1.03 keal/mole and ours 1.09 kcal/mole. The difference is only 6%. 

When the specific heat is known, Debye’s characteristic temperature can be 
computed. This is done with the distribution function at 0°K. As Fig. 4 shows, 
the deviation from the values found by experiment is greater than that found by 
Iona when the temperature is above 10°K. However, he used the distribution func- 
tion at room-temperature, and assumed that it does not change when the temperature 


eA SSeS ee 


1 An attempt has been made to calculate the heights of the peaks by means of the method 
given by van Hove [20] (cf. also Rosenstock [16]). The coefficients are not, however, given as 
analytical expressions in our investigation and the possibilities of estimating an approximate 
expression for the coefficients are limited by the roughness of the network. 
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preliminary calculation of this effect. The contribution of the polarization to the 


coefficients in the equations of motion has been deduced by Lundqvist [13] under — 


the assumption that the polarization can be calculated as if the field were homogene- 
ous, having the same value as at the centre of the ions. Born and Huang [3] have 
considered the general case for an arbitrary lattice. The electronic polarizabilities of 
the ions have here been assumed constant and equal to 1.13 x 10-*4 em’ for K* and 
2.92 x 10-*4 cm’ for Cl-. These values were obtained by Shockley [17]. The contri- 
bution from the polarization to the elastic constants is zero, and so the force constants 
will not change. The only thing we have to do is to add it to the coefficients. We 
cannot, however, regard KCl as a one-atomic lattice, because the polarizabilities for 
K* and Cl- are not the same. 

The polarization seems to have the following effects: the optical branches of the 
frequencies are cut down; the acoustic branches are unchanged for long waves be- 
cause the polarization does not influence the elastic constants, but they are cut 
down for shorter waves. A study of the amplitudes shows that for long waves they 
are equal, but when the waves are shorter they will diverge. Although no definite 
conclucions can yet be drawn, it seems that the frequencies generally decrease and 
the distribution function is pushed to smaller frequencies, and so the specific heat 
will be greater and the Debye temperature lower, in better agreement with the 
experiment. 
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and comments. The work has partly been supported by the Swedish Natural Research Council. 
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